Study of Short-distance Spin and Charge Correlations and Local Density-of-States 
in the CMR regime of the One-Orbital Model for Manganites 
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The metal-insulator transition, and the associated magnetic transition, in the colossal magnetore- 
sistance (CMR) regime of the one-orbital model for manganites is here studied using Monte Carlo 
(MC) techniques. Both cooperative oxygen lattice distortions and a finite superexchange coupling 
among the t2g spins are included in our investigations. Charge and spin correlations are studied. 
In the CMR regime, a strong competition between the ferromagnetic metallic and antiferromag- 
netic charge-ordered insulating states is observed. This competition is shown to be important to 
understand the resistivity peak that appears near the critical temperature. Moreover, it is argued 
that the system is dynamically inhomogeneous, with short-range charge and spin correlations that 
slowly evolve with MC time, producing the glassy characteristics of the CMR state. The local 
density-of-states (LDOS) is also investigated, and a pseudogap (PC) is found to exist in the CMR 
temperature range. The width of the PC in the LDOS is calculated and directly compared with 
recent scanning-tunneling-spectroscopy (STS) experimental results. The agreement between our 
calculation and the experiment suggests that the depletion of the conductance at low bias observed 
experimentally is a reflection on the existence of a PC in the LDOS spectra, as opposed to a hard 
gap. The apparent homogeneity observed via STS techniques could be caused by the slow time 
characteristics of this probe. Faster experimental methods should unveil a rather inhomogeneous 
state in the CMR regime, as already observed in neutron scattering experiments. 

PACS numbers: 75.47.Lx, 75.47.Gk, 75.30.Mb, 75.40.Mg 



I. INTRODUCTION 



The colossal magnetoresistance effect in manganites 
continues attracting considerable interest in both con- 
densed matter physics and material science*^ 2 ^^ Clar- 
ifying the origin of the CMR effect is expected to funda- 
mentally contribute to our understanding of the intrinsic 
complexity observed in a variety of transition metal ox- 
ides. In addition, these compounds are receiving much 
attention for their potential application in new genera- 
tions of spintronics devices^ Manganites present a very 
rich phase diagram with a variety of competing states, 
giving rise to the exotic behavior found in these com- 
pounds. To explain the CMR effect, originally a double 
exchange (DE) model was proposed. 7 However, it was 
shown later that the DE mechanism itself could not ex- 
plain the large changes in the resistivity with increas- 
ing magnetic fields that were observed experimentally. 8,9 
Thus, other couplings, such as the electron-phonon in- 
teractions, must be taken into account^ Subsequent in- 
vestigations have shown that the existence of insulat- 
ing states competing with the DE induced ferromagnetic 
metal, and the concomitant emergence of nanometer- 
scale phase separation, could be the key ingredients to 
understand the CMR effect plil 2 .^ 

After the phase separation idea was proposed, a va- 
riety of theoretical studies have been reported, includ- 
ing calculations using simplified spin models and random 



resistor networks! 14 ' 15 ! 16 More realistic one-orbital and 
two-orbital models, including electron-phonon coupling 
and quenched disorder ) 17 ' 18 i 19 ' 20 i 21 have also been inves- 
tigated. In several of these calculations a CMR behavior 
of the resistivity vs. temperature has been observed. The 
ferromagnetic (FM) metallic phase is commonly accepted 
as one of the competing states necessary for CMR. The 
other competitor should be an antiferromagnetic (AF) 
charge-ordered (CO) insulating state. But for techni- 
cal reasons related with the complexity of the numerical 
simulations, the CMR arising from the competition be- 
tween the FM metallic and the AF/CO insulating states 
has not been discussed in most of the efforts mentioned 
above. However, very recently §en et al3& carefully an- 
alyzed the influence of the superexchange coupling Jaf 
between the localized t2 g electrons on the physical prop- 
erties of both the one- and two-orbital models including 
cooperative oxygen lattice distortions. Two competing 
ground states were discovered at low temperatures: a FM 
metallic phase and an AF/CO insulating phase (for the 
charge and spin arrangement in these states see Fig. [1} . 
These two phases are separated by first-order transitions 
at low temperatures. A thermal transition from either 
phase to a paramagnetic (PM) insulator at high temper- 
ature was also observed. Close to the bi-critical point in 
the phase diagram, the CMR effect in the resistivity was 
found^ 3 - This is a location where enhanced phase com- 
petition is observed. The CMR effect has been shown to 
be deeply connected with the existence of short-distance 



2 



charge correlations above the critical temperature of the 
thermal transition. 22 
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FIG. 1: (color online) The two competing ground states in 
the one-orbital model at electronic density n=0.75, obtained 
by varying the superexchange coupling Jaf, as explained in 
Ref. [22| and in the text. The FM metallic state is in (a) and 
the AF/CO insulating state in (b). The arrows denote the 
classical tz s spins and the open circles are proportional to the 
electronic density. The full circles in (b) simply denote the 
location of the holes. The dashed lines in (b) define "spin 
blocks", that are discussed later in the text. 



In this paper, the recent effort by §en et alM- is further 
expanded, focusing on the analysis of the phase competi- 
tion in the thermal metal-insulator (MI) transition of the 
one-orbital model. We have observed that the develop- 
ment of short-range spin-spin correlations, with charac- 
teristics similar to those in the AF/CO competing state, 
is important to understand the CMR effect. This re- 
sult supplements the previous investigations that focused 
on the charge sector exclusively^ The systems without 
quenched disorder that we study in this work are shown 
to be statistically homogeneous when averaged over very 
long MC times. However, charge localization together 
with robust short-distance correlations of charges and 
spins are observed at short MC times. This leads us to 
argue that a dynamical inhomogeneity occurs in mangan- 
ites in cases where the strength of the quenched disorder 
is weak (namely for materials close to the clean limit). 

A second purpose of our present effort is to use nu- 
merical techniques to calculate the LDOS. Since the one- 
orbital model has been shown to successfully describe 
the behavior of the resistivity in the CMR regime* 2 ^ 
this approach should be able to explain other proper- 



ties in the same temperature range. The LDOS is im- 
portant to contrast theory with STS experiments, which 
are useful to probe local electronic structures in a vari- 
ety of complex materials, possibly providing information 
about phase separation in manganitesi 24 ' 25 In fact, the 
derivative of the experimentally obtained current-voltage 
curve is directly related to the LDOS. Our calculations 
are even more relevant in view of recently reported puz- 
zling STS results. In Refs. [26 ! J27i, a depletion of the con- 
ductance at low bias was observed in the low tempera- 
ture metallic phase of both Lao.77Cao.23Mn03 (LCMO) 
and Lao.35oPro.275Cao.375Mn03 (LPCMO) films, imply- 
ing the existence of a gapped quasi-particle excitation 
spectra in these materials. No coexistence of insulating 
and metallic regions was observed in the experiment, in 
apparent disagreement with existing theories. These re- 
sults are very surprising since a gapped spectrum cannot 
lead to a metallic state. Singh et al. 27 proposed that 
the surface of LPCMO is different from the bulk: at the 
surface the AF/CO state is stable, while in the bulk a 
coexistence of metallic and insulating regions dominates, 
as shown in other experiments. We do not disagree with 
this proposal, and for LPCMO it may well be the solu- 
tion of the puzzle. For LCMO, the same hypothesis may 
work, but here we propose another alternative. In our 
calculations reported below, a PG in the LDOS (as op- 
posed to a gap) is found in the entire temperature range 
where the CMR occurs in the one-orbital model. This PG 
feature is known to be a signature of doped manganites 
and results from mixed-phase tendencies ! 28 ' 29 We have 
investigated the width of the PG and compare it with the 
experimentally determined gap width. We have observed 
a very good agreement between our calculation and the 
experimental results. This provides an alternative way to 
understand the STS data for LCMO: the experimentally 
observed gap could arise from a PG in the LDOS, not 
from a gap. Moreover, within this picture, the appar- 
ent spatial homogeneity observed in the STS experiment 
can be explained as a time-averaged property. The sys- 
tem is actually dynamically inhomogeneous, which could 
be experimentally detected if the time scale of the prob- 
ing technique is comparable to the characteristic time of 
the dynamics of the system. Fast measurements (using, 
e.g., neutron scattering) have already shown evidence of 
"correlated polarons" in LCMO 30 in contradiction with 
LCMO STS experiments. Our calculations show that 
this is to be understood based on the dynamical nature 
of the CMR state: a homogeneous gap detected in STS is 
not in contradiction with other faster measurements, to 
the extent that it is understood that the gap is actually 
a PG in the LDOS. 

This paper is organized as follows: in Sec. II, we briefly 
introduce the one-orbital model and the details of the MC 
technique that we used. In Sec. Ill, we focus on investi- 
gating the competition between the FM metallic and the 
AF/CO insulating states. Spin correlations are shown to 
play an important role in the system, particularly in the 
CMR regime. Moreover, our investigations in this section 
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show that the phase competition in the CMR regime is 
dynamic, at least within the time evolution generated by 
the MC procedure, which is based on local interactions. 
In Sec. IV, a detailed study of the PG in the LDOS 
of the one-orbital model is provided. We compare our 
results with those in the recent STS experiments, and 
show that their observations can be explained within our 
model. Conclusions are in Sec. V. 



II. THEORETICAL MODEL 

In this study, the properties of the one-orbital model 
for manganites are analyzed using numerical techniques. 
Although a two-orbital model with Jahn- Teller couplings 
is more realistic, it is also more difficult to study compu- 
tationally. In fact, using the exact diagonalization proce- 
dure for the fermionic sector, a simulation using two or- 
bitals costs 2 4 =16 times more CPU time than using only 
one orbital. Since it has been extensively shown that the 
one-orbital model captures the physics of phase compe- 
tition in manganitesji 3 - being restricted to one orbital is 
not a drastic approximation. Moreover, the current ef- 
fort builds upon the previous results obtained in Ref. I22I 
where the one-orbital model was also considered. Then, 
for a proper comparison, it is imperative to use the same 
model. The Hamiltonian reads 

H = -t (di ta dj, a + h.c.) - J H ^ d\ )a a a ,pdi,p ■ S, ; 

{ij),a i, a ,{3 

+JAF^Si • Sj - Xt y^u^s - u it s)d\ a di ta (1) 

(ij) z,<5, a 

where di. a and d\ a are the fermionic annihilation and 
creation operators of e g electrons, acting on site i and 
with z-axis spin projection a. {a} denote the Pauli ma- 
trices. Si represents the t2 g spin degrees-of-freedom lo- 
calized at site i. This spin is assumed to be classical, 
a widely used approximation also employed in Ref. I22I 
The first two terms in Eq. [1] correspond to the ordinary 
DE model. The third term introduces AF interactions 
between the localized ^2g spins. Although the coupling 
regulating the strength of this term is considered weak in 
magnitude ^ previous work has shown that its presence 
is crucial to stabilize the AF/CO states that compete 
with the FM ground state resulting from the double- 
exchange mechanism^ 1 . The next term in the Hamilto- 
nian contains the interaction between the mobile elec- 
trons and the cooperative phonons (<5 runs over the spa- 
tial directions). As in Ref. [22| and several other inves- 
tigations, the oxygen lattice displacements Ui t s with re- 
spect to the equilibrium position are considered to be 
classical. These displacements are located at the links 
of a two-dimensional square lattice. The oxygen lattice 
distortions induced by the electron-phonon interaction 



are further balanced by an elastic energy proportional to 
uf s . The last term in the Hamiltonian Eq. [1] contains 
the quenched on-site disorder that simulates the effect of 
chemical doping in real manganites. As in Ref. [22] and 
other publications, this disorder should enter via a ran- 
domly chosen set of values for the on-site energy Aj at 
every site. However, although it has been shown that the 
quenched disorder can enhance substantially the magni- 
tude of the CMR effect, 3 - it has been observed 2 ^ that 
its presence would not change the physical picture qual- 
itatively. If quenched disorder is ignored, the only price 
to pay is the need to tune couplings to be very close to 
the first-order FM- AF/CO transition in order to observe 
CMR physics. 22 Then, for numerical simplicity, here we 
consider only the clean limit CMR, i.e., Aj = and focus 
on the immediate vicinity of the FM- AF/CO transition 
in parameter space. As a consequence, our conclusions 
below regarding the dynamical character of the inhomo- 
geneities will then be of much more relevance to "clean" 
manganites than "dirty" ones. Finally, /1 is the chemical 
potential and rii ta is the number operator. 

In the present calculation, and also for simplicity, the 
limit Jh — ► 00 is considered. This is another of the sev- 
eral well-accepted approximations widely used in previ- 
ous studies of theoretical models for manganites. This 
limit preserves the essential physics of the CMR effect. 
According to this approximation, the spin of the e g elec- 
trons is always parallel to the direction of the spin of the 
electrons. Hence, the Hamiltonian is reduced to 

H = -t / 'X^jjclcj + h.c.) + Jaf / J • Sj 

-Xty^^Ui-s - u it s)c\ci (2) 
i,S 

where c] creates an e g electron with spin parallel to the 
localized t2g spin at site i. The classical localized spin can 
be parametrized as Si = (sin#i cos fa, sin#i sin fa, cos 9i) 
in spherical coordinates. The influence of the infinite 
Hund coupling is reflected in the electron hopping term 
that now carries a Berry phased 

Clij = cos I cos | + sin | sin h^-h) . (3) 

The Hamiltonian Eq.[^]is solved here via a standard com- 
bination of exact diagonalization of the fermionic sector 
and a MC evolution of the classical spin and oxygen lat- 
tice distortion degrees-of-freedom. At each MC step, the 
quadratic fermionic portion in the Hamiltonian is diago- 
nalized numerically, using library subroutines, for a given 
configuration of classical spins and phonons. A new set of 
configurations for the next MC step is then generated fol- 
lowing the Metropolis algorithm. This method has been 
successfully used in previous studies and details can be 
found in Refs. [3llTll and [l3|. Physical quantities can be 
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easily calculated as MC averages. For instance, at each 
MC step, the LDOS is evaluated using the eigenvalues 
and eigenvectors that arise from the diagonalization. 13 
The resistivity p is obtained from the inverse of the aver- 
age conductivity er. In two dimensions (2D), where our 
effort focuses on, it can be shown that a = G, where G 
is the conductance. This conductance is calculated us- 
ing the Kubo formula. 33 In this effort, 2D systems with 
charge density n=0.75 are studied for a proper compari- 
son with the results in Ref. [22| that were obtained at the 
same density and dimensionality. The typical clusters an- 
alyzed here are of size 8x8 and 12 x 12. If not specified, 
simulations start with a random configuration of spins 
and lattice displacements. To overcome the metastabili- 
ties caused by the dynamical inhomogeneity that appears 
in this study (discussed in detail in the following sec- 
tions), very long MC runs were needed. In practice, up 
to 10 5 MC steps for thermalization were used, followed 
by 5 x 10 4 steps for measurements. This was needed on 
the 8x8 lattice to obtain statistically homogeneous re- 
sults. For the same purpose, 10 4 steps for thermalization 
and another 10 4 steps for measurements were used for 
the 12 x 12 lattice. 



III. COMPETITION BETWEEN 
FERROMAGNETIC AND CHARGE-ORDERED 
ANTIFERROMAGNETIC STATES 

It has been shown 22 that the model Eq.[2]has two com- 
peting states at low temperatures, at the n=0.75 density 
considered here. The ground state is a FM metal if Jaf 
is zero or very small. However, upon increasing Jaf an 
AF/CO insulator is reached via a first-order phase tran- 
sition. In this insulating phase, both charges and spins 
are arranged in a non-homogeneous, but regular, pattern. 
This is significantly different from the FM state, where 
both charge and spin are uniformly distributed, namely 
they are the same at every site. At high temperatures, 
the system loses long-range order and it is in a param- 
agnetic (PM) phase. The most striking feature found in 
this model 22 is that close to the phase boundary between 
the FM and AF/CO states, the cooling down of the sys- 
tem from the high-temperature PM phase causes drastic 
changes in the resistivity p. First, p increases with de- 
creasing temperature and develops a prominent peak. It 
then drops rapidly (about 2-3 orders of magnitude) when 
the system enters into the FM state. The resistivity peak 
can be largely suppressed by small magnetic fields, show- 
ing the so-called CMR effect. 22 While it is clear that the 
AF/CO state of the one-orbital model is not exactly as 
found in experiments, which are dominated by CE states 
that also present orbital order, the promising results for 
p vs. temperature show that the key qualitative features 
of the CMR effect do appear in the one-orbital model. 

In this section, we will focus on the MI transition that 
causes the CMR, and show how this CMR effect is con- 
nected to the competition FM vs. AF/CO. 



A. Metal-Insulator Transition and the CMR Effect 

Let us first revisit the thermal transition related to the 
CMR effect that occurs in a narrow parameter regime of 
the one-orbital model. 22 This transition is simultaneously 
insulator to metal and PM to FM, upon cooling, similarly 
as in experiments. This MI transition and the associated 
magnetic transition, for a lattice of size 12 x 12 and at 
electronic density n=0.75, are shown in Fig. [21 The criti- 
cal temperature Tq of the magnetic transition is approxi- 
mately located at the resistivity peak, i.e., Tc ~ Tmi- At 
temperatures lower than Tmi the resistivity falls down 
rapidly, about two orders in magnitude within a narrow 
temperature window T/t sa 0.01. A similar behavior in 
both the resistivity and the uniform magnetization is also 
observed within a relatively narrow range of parameters: 
A - 1.1 - 1.2 and Jaf ~ 0.02 - 0.034. 
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FIG. 2: (color online) Resistivity and uniform magnetization 
upon cooling (black) and heating (red) for Jaf=0.0325 (left 
panel) and Jaf=0.0335 (right panel). The result was obtained 
on a 12 x 12 square lattice using A=1.2, density n=0.75, and 
in the clean limit (A=0). No hysteresis loop is observed in 
both quantities. 

In Fig. [2] the resistivity and uniform magnetization 
phase transitions were studied in the following way. For 
cooling, we started the MC simulation at a high temper- 
ature T 3> Tmi using a random spin and phonon config- 
uration. After MC converging at the initial temperature, 
we started another simulation at a lower temperature by 
using the final configuration from the previous tempera- 
ture as the new initial input. For heating, we proceeded 
in a similar way, but the first simulation at the lowest 
temperature was initialized with a perfect ferromagnetic 
spin configuration (note that the CMR effect occurs very 
close to the FM- AF/CO transition but slightly on the 
side of the FM phase, thus there is no need to start any 
simulation with the AF/CO state). Both the heating and 
cooling procedures were performed "slowly" in temper- 
ature. By this procedure, we can detect possible differ- 
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ent responses of the system upon heating and cooling, 
which is usually associated with the existence of first- 
order transitions. However, based on the results shown 
in Fig. [21 no hysteresis loop has been observed, either in 
p or in the magnetization. This suggests that both the 
Ml and magnetic transitions are of second-order in this 
one-orbital model in the clean limit. Interestingly, our 
result is consistent with a recent analytical work 20 show- 
ing that the transition is of second-order in the absence 
of quenched disorder—. 



B. Short-Range Spin and Charge Correlations 
in the CMR Regime 

The CMR regime, where the resistivity peak develops 
during the MI transition, is close to the bi-critical point 
separating the FM and AF/CO states in the clean-limit 
one-orbital model phase diagram. 22 This indicates that 
at low temperatures, two competing ground states with 
different long-range orders exist. Although by definition 
long-range order cannot survive above a critical temper- 
ature, in this temperature regime it is very common to 
find robust short-range correlations of the order that will 
become stable at low temperatures. This certainly oc- 
curs in the FM phase far from the region of competition 
with other states. However, dramatically different from 
this well-established concept, previous studies 2 ^ observed 
that in the CMR regime, the short-range order that de- 
velops above the Curie temperature is in the charge sec- 
tor. Namely, it is unrelated with the FM state. Moreover, 
the observed short-distance charge correlations are very 
similar to those of the competing AF/CO state. Thus, 
the short-distance physics above Tq near the bi-critical 
point is very different from the properties of the FM state, 
which is the actual ground state of the system. 

The short-distance tendency toward a CO state above 
Tq and Tmi accounts for the appearance of the resistivity 
peak at Tmi. The system develops CO charge correlation 
upon cooling and this proceeds together with an increase 
in the resistance, since the AF/CO state is insulating. In 
this section, we confirm all these previous observations. 
More importantly, here we will also consider the short- 
distance spin correlations above Tq, which were not an- 
alyzed before. We will arrive to the conclusion that both 
short-distance charge and spin correlations are important 
in explaining the resistivity peak. 

An intuitive procedure to understand the short- 
distance spin and charge correlations in the one-orbital 
model is to analyze MC "snapshots" , namely equilibrated 
configurations generated by the MC procedure. Exam- 
ples are presented in Fig. [3] Before analyzing them, and 
to guide the discussion, let us first consider what are the 
expected "building blocks" in the charge-ordered state, 
and then we will investigate if these building blocks are 
observed in the simulations. Since it is anticipated that 
the short-distance behavior above Tc has strong simi- 
larities with the AF/CO competing state, the expected 



building blocks, both for the one- and two-orbital mod- 
els, are shown in Figs. [3] (e) and (f), respectively. In 
the one-orbital model, the study in Ref. [22| has shown 
that at n=0.75 the building block consists of a hole 
carrier surrounded by its four neighboring sites, each 
carrying an electron. The five localized t2 S spins in- 
volved in this block have the same orientation of those 
spins. This block has similarities with a FM polaron, 
but quite contrary to a naive "FM polaron gas" sce- 
nario the AF/CO state is achieved by arranging these 
blocks antiferromagnetically in a regular pattern (see 
also Fig. 1(b)). For the two-orbital model, shown here 
only for completeness, the most likely building block 
should be a small segment of the well-known CE-type 
FM zigzag chains (see Fig. 3(f)). 

It is remarkable that the building blocks of the AF/CO 
state of the one-orbital model can actually be observed 
in the MC snapshots shown in Figs. [3] (a) to (c), for a 
value of Jaf where the ground state is FM. We find that 
the short-range 2 and y/E charge correlations clearly exist 
at short distances, in both the high-temperature insulat- 
ing PM and low-temperature metallic FM phases. In 
the FM state, spins are almost parallel to each other, 
and hole hopping is allowed. This gives the FM state a 
metallic nature. Moreover, it will be shown below that 
the charge is fairly mobile at low temperature, thus ef- 
fectively the building blocks loose their identity and they 
merge into a metallic FM state. Different is the situa- 
tion at larger temperatures. When T > Tmi (shown as 
T//j=0.035 in Fig. [3j, the short-distance tendency to the 
competing AF/CO state is significant. This tendency 
is further confirmed by analyzing the short-range spin- 
spin correlations. On one hand, the spins that are the 
nearest neighbors to a hole site are almost parallel to 
the spin on that site, indicating short-range FM correla- 
tions. On the other hand, spins localized on different hole 
sites prefer to align antiferromagnetically at distances 
\/5, as expected from the discussion in Ref.[22|. Then, the 
MC simulation clearly shows the presence of the building 
blocks described in the previous paragraph (two of them 
are highlighted in Fig. 3(b)). As the temperature cools 
down, these "preformed" building blocks rotate and ori- 
ent themselves into a FM pattern. At higher temperature 
T/i=0.1, both short-range spin and charge correlations 
are much suppressed by thermal fluctuations, but the 
building blocks can still be identified (see Fig. 3(c)). 

By monitoring the above described MC snapshots, a 
simple picture of the temperature evolution of the sys- 
tem near the bi-critical point emerges: upon cooling from 
high temperature, the building blocks develop, and first 
they start arranging in a pattern similar to the AF/CO 
state. This causes the insulating properties above the 
Curie temperature. Upon further cooling, the building 
blocks rotate their relative orientation from AF/CO to 
the FM state, rendering the state metallic. To observe 
this effect more explicitly, we propose an operator O 
which characterizes the spin and charge correlations in 
the AF/CO state, and study the temperature evolution 



6 




T=0.035 

o cD ©A &^>\ ©©© ), © 
0W(!)\(ro GACD 



©/© © ©\CD"V©® 
©©©V ©©© 



j^-©4 ©©©©© 
©0© • ©©© 
©®/©\©©©/© 

©V © 
© \ cp© ©<mm 
©cp i^^gm s 
© t ©©^kk^ 
©©©© T © © 




(e) 



© 

/ N 
/ \ 

<${ $ ^© 



0.00 0.04 0.08 0.12 

T/t 

(f) t 



© 



(1) 



FIG. 3: (color online) (a)-(c) Typical MC snapshots of local 
charge density and spin correlations on a 8 x 8 lattice with 
Jaf = 0.0325, n — 0.75, and A = 1.2 at three different tem- 
peratures. In each snapshot, the circle radius is proportional 
to the local charge density. The 16 holes are represented by 
filled circles. The length of the arrow at each site i is pro- 
portional to the spin-spin correlation (So • Si), where So is 
a reference point located on a hole site. Presenting this cor- 
relation, as opposed to directly showing the classical spins, 
allows an easier visualization since now most spins are in the 
plane of the figure. Holes with y/E and 2 correlations are high- 
lighted, as well as a couple of well-formed building blocks (see 
text) at the two highest temperatures, (d) The temperature 
dependence of the proposed operator O (Eq. |4)| . The result 
of (O) is averaged over 5 x 10 4 MC steps, (e) The building 
block (see text) of the AF/CO state of the one-orbital model 
at n=0.75. (f) Proposed building block of the AF/CO state 
in the two-orbital model at n=0.5, for completeness. 



of O in the system. The operator is 

°= E H^si-sja-mxi 

i.j — hole 



(4) 



where Si = | [Si + J2s ^i+s] is the average of the classical 
spin of the hole site^ and its four neighbors (namely, the 



"spin block" region enclosed by a diamond in Fig. [Ijb) 
or 3(e)). Note that (0)=n^ j B in a perfect AF/CO state 
with nhoic hole carriers. This same operator vanishes in 
either a state with perfect FM order or in a PM state 
with uncorrelated randomly distributed spins. 

The temperature dependences of O for two values of 
Jaf are shown in Fig. [3]^d). At JafA=0.05, the ground 
state does present long-range AF/CO order, and no re- 
sistivity peak was observed. At this Jaf value, (O) in- 
creases monotonically upon cooling and saturates to a 
large finite value when T — > 0. The mean-value of this 
operator illustrates the transition from PM to AF/CO 
states. Interestingly, in the CMR regime, such as for 
Jaf /t=0. 0325 where the FM state is the ground state, 
(O) evolves fairly differently. Upon cooling, it first slowly 
increases following very closely the result at Jaf/£=0.05, 
but then develops a peak at T ss Tmi, and subsequently 
drops fast to a nearly vanishing value at the lowest tem- 
perature investigated. This non-monotonic temperature 
dependence clearly reflects the enhanced tendency to the 
AF/CO state at T > T M i, intuitively observed in the MC 
snapshots. Remarkably, the temperature dependence of 
(O) and the resistivity shows a qualitative similarity. It 
confirms that the resistivity peak in the CMR regime is 
simply a reflection of the short-distance tendency to the 
competing AF/CO state. We note that this idea has been 
proposed in previous studieS ) 21 i 22 but only the charge cor- 
relations were considered as relevant there. It is in this 
work that spin correlations are shown to be equally im- 
portant in understanding the CMR physics that appears 
in the one-orbital model. Charge and spin correlations 
develop together at short distances. 

To better understand the relation between charge and 
spin correlations, the following MC simulation in the 
CMR regime, at Jaf=0.0325 and A=1.2, was carried out. 
The simulation was initialized with a configuration where 
the system has charge correlations similar to the AF/CO 
state. This was achieved by running another simulation 
with J AF =0.05, A=1.2, T/£=0.01, and using the final 
configuration as the initial configuration in the new run. 
A snapshot of this initial configuration is presented in 
Fig. UJa) . During the MC evolution, the oxygen lattice 
distortions degrees-of-freedom were kept frozen to their 
initial values, but the spins were allowed to change (note 
that we use one of the equilibrated MC configurations 
as our "frozen" lattice and there the holes are not regu- 
larly distributed as in a perfect CO state. However, the 
average over many MC configurations at the couplings 
used do give long-range order in the charge correlations) . 
By freezing the phonons, the charge order is stabilized, 
and it provides a static inhomogeneous background that 
is temperature independent. If the spin and charge were 
decoupled, when cooling the system in the CMR regime, 
spins can still experience a PM-to-FM transition, and the 
system may be driven to a metal due to the DE mech- 
anism. However, as shown in Figs. (He) and (d), this 
is not what we obtained from MC simulations. In fact, 
the system is never driven to a FM state even at very low 
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FIG. 4: (color online) MC results obtained by freezing the 
phononic degree-of- freedom (see text for details). Snapshots 
of initial and final configurations are shown in (a) and (b), 
respectively. Also shown are the normalized magnetization 
(in (c)) and resistivity (in (d)). With frozen charges, the DE 
mechanism cannot produce the FM state at low temperatures. 



temperatures. By observing the final configuration of the 
MC simulation at T/t=0.01 in Fig. 0Jb) and comparing 
it with the initial one, we find that spins are effectively 
also frozen if the charges are frozen during the MC evo- 
lution. Thus, the system is locked in the AF/CO state. 
Interestingly, this picture holds even when Jaf=0 (not 
shown). Then, clearly spin flipping processes are greatly 
suppressed by freezing the charge and lattice distortions. 
This suggests that: (i) the spin and charge degrees are 
strongly coupled in the CMR regime; (ii) CMR cannot 
take place in a quenched inhomogeneous charge back- 
ground. The CMR is related to a more complex effect: 
upon cooling across Tq, the charge order partially melts 
allowing for a metal to develop. 



C. Dynamical Properties of Charge and Spin 
Correlations in Monte Carlo Time 

In the one-orbital model without quenched disorder be- 
ing analyzed here, the translation invariance symmetry 



cannot be broken unless a state with long-range AF/CO 
order is stabilized, and this occurs only at low temper- 
atures and sufficiently large Jaf- Then, in the CMR 
regime where the resistivity peak appears, the MC-time- 
averaged charge density must be the same at every site. 
However, in the MC snapshots previously discussed, a 
short-distance charge localization resembling the AF/CO 
state leads to a non- uniform distribution of charge (rij). 
Then, as a state without long-range CO but contain- 
ing strong distortions resembling the AF/CO state, the 
charge localization in the CMR regime can only be under- 
stood in a dynamical context. In other words, even after 
the MC thermalization that removes the dependence on 
the initial configuration, in the CMR regime the system 
is still inhomogeneous at any instant during the simula- 
tion, but it slowly becomes homogeneous when averaged 
over very long MC times. It is, thus, interesting to in- 
vestigate this dynamics. Note that, in principle, it is 
not obvious that real-time and MC-time dynamics have 
any common trends. However, both are based on local 
events (local interactions in the real system to evolve in 
real time; local updates in the MC evolution used here). 
Thus, the expectation is that long-lived metastabilities 
and tendencies toward glassy behavior and a complex 
landscape revealed via MC-time studies will have an ana- 
log in the real system. To the extent that the analysis in 
MC time is considered only qualitatively, our expectation 
is that crude dynamical trends for the real system can be 
discussed based on a local-update MC-time evolution. 36 
With these caveats, here several characteristic MC times, 
revealing dynamical correlations of the spin and charge 
degrees-of-freedom in the system, will be numerically an- 
alyzed. It will be argued that glassy properties emerge 
at the CMR regime, and we speculate (but do not prove 
rigorously) that similar tendencies will occur in real man- 
ganites. A variety of experiments showing glassy tenden- 
cies indicate that this assumption is realistic. 37 

Consider first the MC-time evolution of the conduc- 
tance (inset of Fig. [5]). In the present type of studies 
based on calculating the quantum-mechanical transmis- 
sion to obtain conductances, it is well-known that the 
conductance is quantized, namely it varies among dis- 
crete values due to the presence of an integer number 
of channels of transmission, from one side to the other 
of the cluster investigated. The observed "jumps" be- 
tween two quantized values reflects on changes that oc- 
curred in the charge and spin correlations. Thus, the 
average MC time t& that the conductance remains at a 
certain quantized value provides information about the 
system MC dynamics. From Fig. [51 it is interesting to 
observe that tc has a very clear peak precisely at the 
metal-insulator transition temperature. In other words, 
as the system is cooled down, concomitant with the de- 
velopment of short-range charge and spin correlations, a 
variety of metastable states develop, as in a glass. The 
system tends to be trapped in those states and evolving 
to other states is difficult. Eventually at very low tem- 
peratures, deep into the FM state, the metastabilities 
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FIG. 5: (color online) Characteristic times that illustrate the 
spin and charge MC dynamics in the one-orbital model are 
here shown as a function of temperature, together with the 
resistivity, to is the average MC time that the conductance 
G remains constant at one of its quantized values. r c is the 
autocorrelation time of the local charge density. r s is the 
autocorrelation time of the uniform spin structure factor. The 
inset contains typical measurements of the conductance as a 
function of MC time, illustrating the discrete values it takes. 



vanish and the system rapidly evolves in MC time. It is 
remarkable that to qualitatively follows the shape of the 
resistivity p. This suggests that the charge localization 
close to Tmi has glassy characteristics^ and if the free 
energy were available in the CMR regime, it would show 
a complex landscape of peaks and valleys. The long life- 
times near the CMR resistivity peak show the practical 
difficulty in achieving true translation invariance: for this 
to occur, the system must visit dynamically all the many 
competing arrangements of charge and spin to render the 
average uniform. 

We can also study an autocorrelation time r c related 
with the charge. For this purpose, let us define the au- 
tocorrelation function A(t) of the local charge density as 



A(t) = 



1 V T ° \<n t+t> 

NT £-ii=\ £-it' = \\."-i 



(•5) 



where rij = tjt J2t=i n \i ano - is the total number of MC 
steps used in the measurements (50,000 in this case). It 
was observed that A(t) decays exponentially as A(t) ~ 
e _ */ T ° , where t c is the autocorrelation time of the local 
charge density. The temperature dependence of r c is also 
shown in Fig. It has a similar shape as tQ > indicating 
that charge has a tendency toward localization in long- 
lived states near the CMR peak. However, the peak in r c 
is shifted to temperatures lower than the metal-insulator 
transition temperature. Thus, charge localization is still 
quite stable even below Tmi- This picture changes at 
very low temperature, where r c increases again, and the 



system stabilizes a metallic state with uniform charge 
distribution. This increase of r c at low temperature is 
unrelated to charge localization, but it reflects on the 
stable uniformity of the FM state charge distribution. 

As discussed above, spin dynamics is also important in 
understanding the resistivity peak in the CMR regime. 
Similar to charges, an autocorrelation time can be defined 
for the spins. This autocorrelation function reads 



A(t) = 



TlLi[s t+t ' -s][5*'-s] 



(6) 



where S = jj J2i j &i ' &j 1S the uniform spin structure 
factor, and S = -jr Y^tli T ne autocorrelation time 



t s is then defined through A(t) ~ e 



-t/T„ 



As shown in 



Fig. [5l it increases quickly at T < T MI indicating the 
stabilization of a FM state. 

We found that very short-range FM correlations ex- 
ist even in the AF/CO state. This is clear by observ- 
ing a typical MC snapshot of spin correlations. In a 
perfect AF/CO state, each hole is surrounded by four 
neighbors that have their classical spins parallel to the 
classical spin on that hole site. As previously discussed, 
such spin blocks are then antiferromagnetically arranged 
in the AF/CO state (see Fig.[IJb)). However, if all these 
block spins align, a FM state is reached. It is then in- 
teresting to ask how the spins rotate when experiencing 
a transition from the short-distance AF/CO state in the 
CMR regime to the FM state stable at lower temper- 
atures. Do they rotate coherently, behaving as blocks 
of spins, or the spins rotate more individually one at a 
time? To clarify this matter, let us study the MC evo- 
lution of the average block spin |S| = J ^ J2i • Let 
us also study the average angle between the spin at a 
hole site and its four nearest-neighboring spins, defined 
as cosO = (1 _» )£ 2 J2i=hoies s i ' s i+*- Tne results are 
shown in Fig. [SJ The MC simulation was performed 
at low temperature where a FM ground state is finally 
reached. However, it was started with an initial config- 
uration obtained in an almost perfect AF/CO stated 
The initial and final configurations have very similar val- 
ues for |S|. If the spins rotate incoherently, |S| should be 
first reduced substantially and then grow back when a 
large portion of spins are parallel (to establish the long- 
range FM order) , and cos O should be randomly centered 
around 0. If, on the other hand, spins rotate coherently, 
both |S| and cosO should remain nearly constant and 
close to 1 over the MC time evolution. As seen from 
Fig. Efa), both quantities converge to their equilibrium 
values very fast, even before the thermal equilibrium of 
the FM state is established. It is interesting that both 
quantities do not change too much during the MC time 
evolution, which indicates that spins rotate coherently in 
blocks, at least as a first approximation. 

These results would suggest the following picture: al- 
though not perfectly, the spins in a block mainly rotate 
coherently, such that they keep parallel to one another 
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FIG. 6: (color online) (a) MC time evolution of \S\, cos0, 
and the normalized magnetization m. The MC simulation 
is initialized by a configuration with short-distance AF/CO 
order (similar to that shown in Fig.fjjb)), but the simulation 
occurs at very low temperature, where a FM ground state 
is stable (similar to that shown in Fig. HJa)). (b) MC-time 
evolution of \S\ and cos0, at shorter MC-time scales than in 
(a). Three snapshots of a typical block of spins in the early 
stages of MC evolution are also shown. 



during the MC evolution. This picture is consistent with 
the charge localization scenario for the AF/CO state, and 
the MI transition in the CMR regime can be understood 
in the following way. At T > Tmi, holes are localized 
in these local FM spin blocks, and different spin orienta- 
tions between adjacent blocks suppress inter-block charge 
transfer, giving the PM phase an insulating nature. How- 
ever, holes are not frozen at their initial sites. As they 
travel from site to site, they coherently rotate spins on 
their neighboring sites in the block. At T<T C « Tmi, 
long-range spin correlations establish. Blocks are aligned 



in parallel, charges are no longer localized in blocks but 
can freely move all over the lattice, developing metal- 
lic properties. Note the consistency between the above 
described picture and the "correlated polaron" picture 
often mentioned in the experimental literature^ by re- 
ferring to the local FM blocks around holes as "small po- 
larons" . These polaronic entities are not independent but 
they have tendencies to very particular short-distance ar- 
rangements in the CMR regime. This agreement between 
theory and experiment also supports the view that at 
least qualitatively the study of MC dynamics provides 
useful information regarding the actual real time dynam- 
ics of manganite compounds. 



IV. THE LOCAL DENSITY-OF-STATES IN THE 
ONE-ORBITAL MODEL 

Motivated by recent STS experiments^^ we have 
systematically studied the behavior of the LDOS in the 
one-orbital model for manganites, with focus on the CMR 
regime. For this analysis, we have worked on an 8 x 8 lat- 
tice with model parameters Jaf=0.0325 and A=1.2. Be- 
fore proceeding with the local DOS analysis, note that 
the site averaged LDOS is just the DOS N(w) of the 
system. Unlike what the experimental results have sug- 
gested recently,— we do not observe any hard gapped 
feature in N(w), at any temperature. However, as shown 
in Fig. [7J in the CMR regime of the one-orbital model, 
our results for the DOS present a clear and fairly deep 
PG centered at the chemical potential, both below and 
above the critical temperatures. Starting from the high- 
temperature insulating phase, the DOS at the chem- 
ical potential first decreases upon cooling the system, 
it reaches a minimum at approximately Tmi, and then 
grows up again. The temperature dependence of the in- 
verse of the DOS, N" 1 ^ = //), qualitatively follows the 
shape of the resistivity. An enhancement in iV _1 (w = //) 
reveals a clear tendency to an insulating state, result 
compatible with all the previous analysis of transport 
and spin/charge short-distance correlations. Once spa- 
tial fluctuations are included, most features observed in 
the DOS apply to LDOS as well. 

Although a hard gap in the LDOS was not observed in 
our model even in the insulating phase, we found sev- 
eral similarities between properties of the PG regime 
and those of a normalized conductance curve reported 
in experiments^ The obvious one is that the PG exists 
in the full CMR temperature range we studied, both in 
the FM metallic state and in the PM insulating phase, 
while in experiments a hard gap is found in both regimes 
as well. Moreover, in the experiment, the width of the 
LDOS gap near the chemical potential \l was determined. 
Its temperature dependence was associated with the de- 
velopment of a polaron binding energy. For further com- 
parison with experimental results, we have studied the 
width of the PG in our LDOS. The local width AE t 
is determined by taking the half-distance between the 
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FIG. 7: (color online) Left: Density of states N(u)) at various 
temperatures, in the CMR regime. A pseudogap develops 
in both the PM insulating and FM metallic states. Right: 
temperature evolution of N^ 1 ^ = fx) and the resistivity p. 



two peaks that are the closest to the chemical potential, 
at each lattice site i. The site-averaged width AT, its 
variance cab, and the distribution of AT.;, denoted by 
P(AEi), are then calculated. 
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FIG. 8: (color online) Temperature dependence of the site- 
averaged width AE of the pseudogap in the LDOS. Also 
shown is the resistivity (black diamonds) at the same model 
parameters. Inset: AE rescaled by the critical temperature 
of the MI transition Tmi ■ Also shown as black squares are the 
experimentally determined "polaronic gaps" (extracted from 
Ref.[13) rescaled by Tmi. 



In Fig. [51 the site-averaged width AE is shown to- 
gether with the resistivity. AE has a small value at 
low temperatures in the FM metallic phase, as expected. 



However, this quantity increases upon raising the tem- 
perature. Close to T ~ Tmi, AT develops a cusp and 
it flattens into a plateau at about where the resistivity 
shows a peak. Actually, at exactly T w Tmi it reaches a 
local minimum, namely a little dip on the plateau, but it 
is difficult to judge if such a small effect can survive the 
bulk limit. Upon further heating, AE slowly increases 
again. To directly compare our results with experimen- 
tal data, we rescaled AT and the temperature T with 
respect to Tmi- The rescaled dimensionless AT together 
with experimental data (rescaled as well) are shown in 
the inset of Fig. [SJ The quantitative agreement is very 
good in almost the full temperature range. Nevertheless, 
there is a difference between experiments and our calcu- 
lation close to T w Tmi- Experimentally a notorious dip 
was observed, but this feature appears to be very weak 
in our simulations, as already discussed. However, the 
overall behavior is at least qualitatively the same both in 
theory and experiments. 
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FIG. 9: (color online) Statistics of the width of the LDOS 
pseudogap AEi. (a) A general Gaussian probability distribu- 
tion of AEi at various temperatures in both the FM metallic 
and PM insulating phases is observed, when MC-averaged 
over long times, (b) The standard deviations of AEi. The 
system behaves more inhomogeneously in shorter MC-time 
scales. 



As a further comparison, we present the distribution 
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of LDOS widths, P(AE), in the PM insulating phase, as 
well as in the FM metallic phase. As shown in Fig. [Ha), 
we have observed that in both phases the distribution 
can be fitted by a Gaussian, in full agreement with the 
experimental findings^ At first glance, a single-peak 
Gaussian distribution implies that the system is homo- 
geneous, with no competing metallic and insulating do- 
mains. However, we have to emphasize that the distribu- 
tion in our model is taken from very long MC simulations, 
similarly as the STS experiment which is "slow" in typi- 
cal time scales of relevance for electronic systems. Thus, 
the Gaussian distribution of widths indicates that the 
system is statistically homogeneous when time averaged. 
In other words, the site-averaged LDOS PG width, AE, 
equals to the PG width of the DOS in a very long MC 
run. But it does not exclude the possibility of having 
a dynamical inhomogeneity, such as the charge localiza- 
tion and formation of local FM spin blocks discussed in 
the previous sections. Since the distribution of the PG 
width is Gaussian, its variance ctae is then a measure of 
the inhomogeneity in the system. This quantity is illus- 
trated in Fig.[HJb): ctae strongly depends on the number 
of MC steps in the simulation, especially close to Tmi in 
the insulating phase, where the dynamical charge local- 
ization and correlations are enhanced. If the number of 
MC steps for the MC measurements is comparable to the 
characteristic time of the system dynamics, such as the 
autocorrelation times t c and r s , a different distribution 
of LDOS PG widths may appear. 

Following this idea, we compare the distribution of 
AEi between two MC simulations. These simulations are 
performed for the same couplings and other parameters, 
but one is done with 5 x 10 4 MC steps, while the other 
one has only 100 MC steps (both after 10 5 MC steps for 
thermalization). As shown in Fig. [TUT a). the distribu- 
tion is now quite different: a Gaussian is observed for 
the long run, but a double-peak distribution, which can 
be fitted by a superposition of two Gaussian's, applies 
to the short run. The real space images in Fig. fTUlb) 
and (c) further confirm that the system looks more inho- 
mogeneous at very short time scales. For the same rea- 
son, if one were able to perform the STS experiment in a 
much shorter time scale than currently possible, namely 
with times comparable to characteristic relaxation times 
in the system, then a different gap distribution revealing 
the intrinsic dynamical inhomogeneity would be observed 
in the experiment. 

It is important to further discuss this surprising agree- 
ment between our calculation and the STS experiments. 
While it is easy to imagine a gapped spectrum in an in- 
sulator, finding a gap in the metallic phase seems very 
unusual. The explanation proposed in Ref. [2?] is clearly 
a possibility, namely the surface behaves differently than 
the bulk. However, if the concept of a hard gap is re- 
placed by the pseudogap, then another possible scenario 
appears. A PG in the metal renders the system a "bad 
metal" , as observed experimentally since the values of the 
resistivity are high in the low temperature phase. Find- 
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FIG. 10: (color online) (a) Distribution of the width of the 
LDOS pseudogap AEi at different MC time scales. Real space 
images of AEi measured under 50, 000 and 100 MC steps at 
T/t — 0.02 are also provided in (b) and (c), respectively. The 
radius of each circle is proportional to the local AEi value at 
that site. Open circles corresponds to AEi/t ^ 0.2, whereas 
solid (red) ones corresponds to AEi/t < 0.2. 



ing charge approximately localized and correlated in the 
FM metallic state near the MI transition is not too sur- 
prising, according to the physical picture and results we 
have discussed in previous sections. It is then natural to 
explain the STS gapped feature of the conductance spec- 
trum observed experimentally as a result of a fairly deep 
PG, which cannot be distinguished experimentally from 
a hard gap. In retrospect this is fairly obvious: a hard 
gap cannot exist in a metal. But for a PG coexisting 
with bad metallic properties, there is no contradiction. 



V. CONCLUSION 

In conclusion, we have investigated the MI transition 
in the CMR regime, and the associated PM to FM tran- 
sition, using the one-orbital model with cooperative oxy- 
gen lattice distortions and superexchange coupling Jaf, 
at electronic density n — 0.75. Both the resistivity and 
magnetization are continuous as a function of tempera- 
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ture, indicating a second order phase transition. How- 
ever, strong tendencies to develop spin and charge short- 
range order were observed. It has been shown that the 
tendency toward the AF/CO state at short-distances ac- 
counts for the appearance of the sharp resistivity peak 
at the critical Curie temperature. Robust charge local- 
ization and correlations exist in both the insulating and 
metallic phases, and are enhanced in the temperature 
range where a resistivity peak exists. Charge localiza- 
tion together with correlations in the spin and charge 
degrees-of-freedom are shown to be dynamical, at least 
with respect to the MC time in the simulations, which 
can be characterized by autocorrelation times t c and r s . 
The model studied here presents properties in the CMR 
regime that resemble a glassy system. 

We have also studied the LDOS of the one-orbital 
model in the CMR regime. A PG is observed at all tem- 
peratures. The width of this PG, together with the vari- 
ance and distribution of the width, are calculated and 
compared with results from a recent STS experiment. 
To our surprise, given the crudeness of the calculations, 
theory and experiments agree not only qualitatively, but 
in some cases even quantitatively. This study provides 
an explanation to the puzzling STS experimental results, 
suggesting that the gapped features in the normalized 
conductance spectrum observed in the STS experiment 



actually originates from a pseudogap. We also point out 
that the observed homogeneity in STS investigations may 
reflect on the time-averaged properties of the system. 
The system is actually dynamically inhomogeneous, at 
least in our Monte Carlo time evolution. Our predic- 
tion is that experimental techniques that are fast, namely 
with time scales comparable to the system's intrinsic dy- 
namics in the CMR regime, should see indications of the 
widely expected nanoscale phase separation. However, 
slow techniques, such as STS, should observe an homo- 
geneous state in the same regime. 
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